QCD sum rules for J/ip in the nuclear 
medium: calculation of the Wilson coefficients 
of gluon operators up to dimension 6 



Sungsik Kim and Su Houng Lee 

Department of Physics and Institute of Physics and Applied Physics, Yonsei 
University, Seoul 120-749, Korea 



Abstract 

We calculate the Wilson coefficients of all dimension 6 gluon operators with non 
zero spin in the correlation function between two heavy vector currents. For the 
twist-4 part, we first identify the three independent gluon operators, and then pro- 
ceed with the calculation of the Wilson coefficients using the fixed point gauge. 
Together with the previous calculation of the Wilson coefficients for the dimension 
6 scalar gluon operators by Nikolaev and Radyushkin, our result completes the list 
of all the Wilson coefficients of dimension 6 gluon operators in the correlation func- 
tion between heavy vector currents. We apply our results to investigate the mass 
of J ftp in nuclear matter using QCD sum rules. Using an upper bound estimate on 
the matrix elements of the dimension 6 gluon operators to linear order in density, 
we find that the density dependent contribution from dimension 6 operators is less 
than 40% of the dimension 4 operators with opposite sign. The final result gives 
about —4 MeV mass shift for the charmonium at rest in nuclear matter. 
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1 Introduction 



Identifying higher twist operators and calculating their corresponding Wil- 
son coefficients are very important in several aspects. First, these provide a 
systematic building blocks to analyze the available data from deep inelastic 
scattering(DIS) at lower Q 2 region. Second, these contributions are essential 
for a realistic generalization of QCD sum rule methods to finite baryon density. 

The twist-4 operators were classified and its anomalous dimensions for some 
of the operators were calculated by S. Gottliebfl]. However, the operators 
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were over-determined and the independent set of twist-4 operators appearing 
in the DIS were first identified by Jaffe and Soldate[2], where the twist-4 
operators are of four quark type and quark gluon mixed type. Among the uses 
of these results, the available DIS data were analyzed to determine the nucleon 
matrix elements of the twist-4 operators[3,4]. These estimates were then used 
in generalizing the QCD sum rule approaches for light vector mesons to finite 
density [5, 6]. 

In the correlation function of two heavy vector currents, only gluon operators 
contribute in the operator product expansion(OPE). This is so because in 
the heavy quark system, all the heavy quark condensates are generated via 
gluonic condensates [7-9]. In dimension 6, there are scalar operators, twist-2 
and twist-4 operators. For the scalar gluonic operators at dimension 6, there 
are two independent operators. In ref.[10], the two were identified and the 
corresponding Wilson coefficient were calculated . For twist-2 gluon operator, 
the calculation for the leading order(LO) Wilson coefficient is simple and its 
matrix element is just the second moment of the gluon structure function. 

The twist-4 dimension 6 gluon operators are more involved. In this work, 
we have identified the three independent local gluon operators and calculated 
their corresponding LO Wilson coefficients in the correlation function between 
two heavy vector currents. This result is new and complimentary to a previous 
work[ll] on gluon twist-4 operators, where they start from certain diagrams 
and identify the three independent twist-4 gluon structure functions. Together 
with the previous calculation of the Wilson coefficients for the dimension 6 
scalar gluon operators by Nikolaev and Radyushkin[10], our result completes 
the list of all the Wilson coefficients of dimension 6 gluon operators in the 
correlation function between heavy vector currents. As an application, we will 
use our result in QCD sum rule approach to investigate the property of J/ip 
in nuclear matter. 

This is particularly interesting because the on-going discussion of J/ip sup- 
pression in RHIC as a possible signal for quark gluon plasma[12], inevitably 
requires a detailed knowledge of the changes of J/ip properties in "normal" 
nuclear matter [13]. Furthermore, the large charm quark mass m c ^> Aqcd 
provides a natural renormalization point for which a perturbative QCD ex- 
pansion is partly possible. In fact, studies have shown that the multi-gluon 
exchange between a cc pair and nucleons might induce a bound cc state with 
even light nuclei [14-20]. In such analysis, the low energy multi-gluon potential 
was modeled either from the effective theory obtained in the infinitely large 
rriQ limit [16-18] or from extrapolating the high energy scattering via pomeron 
exchange to lower energy [14]. Although both approaches, gave similar binding 
for the J/ip in nuclear matter, it is not clear how reliable these results are 
unless one systematically calculates the corrections. 
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In order to confirm this findings in an alternative but a systematic approach 
and to establish a basis for further studies, we have previously applied the 
QCD sum rules [7, 8] to heavy quark system in nuclear medium and calculated 
the mass of J/ip and rj c in nuclear medium [2 1,2 2]. This was the generaliza- 
tion of the sum rule method for the light vector mesons in medium[5] to the 
heavy quark system. It was found that the mass of J/ip (i] c ) would reduce 
by 7 MeV (5 MeV), which is indeed consistent with previous results based 
on completely different methods[15,16,19,20]. However, it was not possible to 
reliably estimate the uncertainties of the result, because the contribution from 
the operator product expansion was truncated at the leading dimension 4 op- 
erators. To overcome the limitations , we will here make use of our calculation 
to include the complete dimension 6 contributions. Unfortunately, at present, 
there is no data to identify the nucleon expectation value of any of the di- 
mension 6 operators(scalar, twist-2, twist-4) nor is there any lattice result. 
Nevertheless, in this article, we will estimate the nucleon matrix elements of 
these operators and then use QCD sum rule approach to study the reliability 
of our previous result on the mass shift of J/ip in nuclear matter. If in the 
future, the matrix elements are better known, we will be able to study the 
non trivial momentum dependence of the J/ip in nuclear medium, as has been 
done for the light vector mesons[6]. The momentum dependence is especially 
interesting because, there are inelastic channels opening when the charmo- 
nium system is moving with respect to the medium. For example, when the 
charmonium is moving with sufficient velocity, it will have enough energy to 
interact with a nucleon to produce a D meson and a charmed nucleon. This 
effect would be very important in relation to J/ip suppression in RHIC. 

In section 2 we characterize all gluon operator up to dimension 6 and obtain 
identities to be used to reduce the operators to an independent set. In sec- 
tion 3, we calculate the Wilson coefficients for the independent set of gluon 
operators and show current conservation. In section 4, we calculate the mo- 
ments and perform a moment sum rules analysis to calculate the J/ip mass 
in nuclear matter. We conclude with some discussions. The appendix includes 
some detailed calculation of the Wilson coefficients. 



2 Operators 

In the operator product expansion of heavy quark system, only gluonic con- 
densates are relevant. This is so because all the heavy quark condensates can 
be related to the gluon condensates via heavy quark expansion [7-9]. This is 
also true in nuclear medium since there are no valence charm quarks to leading 
order in density and any interaction with the medium is gluonic. Let us start 
by categorizing gluonic operators up to dimension 6, which does not vanish in 
nuclear matter. 
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For dimension 4 operators, the scalar and twist-2 gluon operators contribute 
[21,22], 



9 2 G a ^ u G a ^ u , g 2 G a fJLa G a ua . (1) 



For dimension 6 operators, one can think of generating a number of gluon 
operators constructed from three gluon fields G a or two gluon fields with 
two covariant derivatives. However, for scalar operators, there are only two 
independent scalar operators [10]. They are, 



3 eabc/~ta s~ib s~ic 2/-ia s~ia /<->\ 

y J & imx va;vfj,' \ z ) 

The second operator can also be written in terms of four quark operator using 
the equation of motion 



G%. v = gqi^—q = gj;, (3) 



As for the spin 2 operators in dimension 6, which are also called dimension 
6 twist 4 (twist=dimension-spin) operator, we can first categorize possible 
operators as follows. First, there is again one operator with three gluon field 
strength tensor. Then, assuming the free symmetric and traceless indices to 
be n and u, depending on whether or not the free index goes into the covariant 
derivative, there are 6 more operators possible. Starting with the three gluon 
operator, the 7 are, 



3 rabc/^ia /~ib /~ic 3 £/^3 

g j ^vkWa^ka = g i^fiu 

g ^kX^k\;ij,v 

g u /j/t u 'v\,k\i g f/,n v^wi g ^^^faiak 

g 2 G° lLK G a KX . Xv , g 2 G a KX G° iiK . Xv . (4) 



However, they are not independent. Using the identities in appendix A, one 
can show that there are three independent spin-2 operators at dimension 6. 
The set we will use are, 



g 2 GkxG^x^v, g 2 g® k gi x . Xk , g 2 g^ k g^, x . Xu . (5) 
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It is interesting to note that in reference [11] starting from certain diagrams 
with two, three and four gluon exchange in the t-channel, they were able 
to derive three twist-4 gluon distribution amplitudes, from which one can 
calculate the nucleon matrix element of the three independent operators. 

As for the spin 4 dimension 6 operator, it is just the twist- 2 gluon operator. 



3 Polarization (OPE) 

Having established the independent gluon operators in dimension 6, we will 
calculate their LO Wilson coefficients in the correlation function between two 
vector currents made of heavy quarks, j p = h'j^h. 



where {-) p represents the expectation value at finite nuclear density p. The 
fourier transforms are defined by 



S(x, 0) is the heavy quark propagator in the presence of external gluon field 
[23]. To calculate the Wilson coefficients, we obtain the quark propagator in 
the presence of external gauge fields in the Fock-Schwinger gauge [10,23], 



In appendix C we list the momentum space representation of the quark prop- 
agator multiplying the gauge invariant gluon operators. We will use these 
quark propagators in eq.(6) and extract the Wilson coefficients by collecting 
the appropriate tensor and gluon structures. 

In general the polarization tensor in eq.(6) will have two invariant functions. 
They can be divided into the longitudinal and transverse part of the external 
three momentum, which are the following when the nuclear matter is at rest. 




(6) 




(7) 



x,A»(x) = 0. 
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where q = (cu, k) is the external momentum. 

In the vacuum or in the limit when k — > they become the same so that there 
is only one invariant function 



n L (^ 2 ) = n T (c 2 ) = ^<nv = n(c 2 ). (10) 

So in this work, we will construct the sum rule for n(c<J 2 ) at k — > and nuclear 
matter at rest. However, in the calculation of the OPE, we will start from the 
general expression of eq.(6) at nonzero value of k and calculate each tensor 
structure separately. This way, current conservation will be a nontrivial check 
of our calculation and the generalization to k ^ will be straightforward [24]. 

In the following subsections, we will summarize the OPE for operators of 
dimension 4 and 6. The new results of the present work are the OPE for 
dimension 6 spin 2 and dimension 6 spin 4 operators. 



3.1 scalar contributions 



Here we summarize the OPE of scalar operators of dimension 4 and dimension 
6 used in the vacuum sum rule[10]. 



np lar (g) = (q,q„ - g»„q 2 ) [C pert - + (g G a al3 G a al3 ) 

+C° G s(9 3 f abc G a a AG c Xa ) + C^(g 4 f K f K )}, (11) 

where 



^ = ^4¥ [ - 1 + 3J2 " 2J31 ' 



CO 



C°. 



G3 727I- 2 (Q 2 ) 3 
1 



' 36tt 2 (Q 2 ) 3 

and 



2 1 31 43 12 

y + 4J 2 J3 H J4 J5 

.15 10 3 5 5 

41 3 ,2 1 . 7 4 7 26 7 8 T 

is" - 5 V + ( 3 + 3 V)Jl - J2 ~ 9 J3 " W J4 + 5 J5 



(12) 



^(y) = / r - ,f : wv , (13) 
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with y = Q 2 /m 2 and m being the heavy quark mass. This will give the fol- 
lowing contribution to II (u 2 ) defined in eq.(10). 



u(uJ 2 ) = C per t . + Cq 2 {g 2 G a a pG a aj3 ) + C G3 {g s f^G a ap G^G c Xa ) 

+C^(g 4 f K f K ). (14) 

For later convenience, we will use g 3 f abc G a ^G b uX G c X)1 = \g 2 G a ^ v G a ^. aa + g*j£j° 
and rewrite the OPE as follows, 



U(u 2 ) = C pert - + C G 2(^-G a tll/ G a lxv ) + C GD 2 G (^Gl u G a ^. KK ) 
2 n 

+C P --(-G a aK G a aX . M ), (15) 

6 71 



where 



^ = ^[-1 + 34-2.7,], 

C GD 2 G = 108Q | g2)3 [4 -3y + 120 J 2 - 310J 3 + 258J 4 - 72 J 5 ], 

C P = 108Q | g2)3 [(4 - 3y + 120J 2 - 3IOJ3 + 258J 4 - 72J 5 ) 

-180 + 120y - (120 + QOy) J x - 300J 2 + 1320J 3 - 720 J 4 ]. (16) 
Note Q 2 = —uj 2 in our kinematical limit. 



3.2 dimension 4 spin 2 operator 



This part has been calculated recently [21,22]. Here we summarize the results. 



rt 2 (g) = ^ 



J 2 



Q 2 



r2 \ 
-pp.) 



+g^ 



q p q* 



q 2 -J pa 



Q 4 



pa 



+ Jpl) 



(17) 



where 



pv 



CzG^Gl v ) 

71 



J, 



{—G'LG'L 

7T 



crp 



(18) 
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3.3 contribution from dimension 6 and spin 4 



The diagrams describing interactions with the gluonic field in dimension 6 are 
shown in Fig 1. 

The calculation for this part is straightforward. We substitute into eq.(6) parts 
of the quark propagator containing D's and G"s (summarized in appendix 
C) so that when these from the two propagators are combined, yield terms 
proportional to GDDG. In Fig.l, these come from Fig. 1(c) to Fig. 1(g). Then 
we extract the traceless and symmetric spin-4 part of the operator GDDG. 
The final result for the dimension 6 spin 4 part of the operators yields, 



n^ 4 (g) 



2\3 



^t,Xpu + Q2 i^P^KXpv + 'ip^nXpp) 



(Q 

QpQa j4 _|_ QpQvQpQv I r4 
^ U Q2 KXpo + Q4 V 1 ! 



+ T 4 ) 

nXpa ' nXpa ) 



(19) 



where 



I* = 



pupa 



266 



45 

x/ — C a C a \ 

> pK , <TK\IXV I 



20 22 



7T 



138 916 16 

Ji — t/a ~\~ — Ja 

5 45 3 



J 4 

pi/pa 



362 / 22 \ 94 T 44 T 16 T 32 T 

4 H w Ji J 2 J 3 H J 4 J 5 

45 V 15 y / 15 45 3 15 



,a 



X { — ^pK^aK; 



71 



,;pvl 



(20) 



3.4 contributions from dimension 6 and spin 2 



The calculation for this part goes in a similar way as in the previous subsection. 
However, this time all the graphs in Fig.l contribute. Moreover, we have to 
include all possible combination of G and D's that makes up dimension 6. We 
then extract the spin 2 part of each operator and use the identities given in 
the appendix to reduce the operators into the independent three in eq.(5). An 
example needed to extract the spin 2 part from a general operator is given in 
appendix D. Finally we find, 



KM 



(Q 2 y 



l lv + qiiw^pv + QpQuI 2 ^) 
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Fig. 1. 1-loop diagrams describing interactions with the gluonic field giving dimen- 
sion 6 operators: The dashed line represents an external gluon field G and the dashed 
line with n arrows represents an external gluon field with n derivatives D n G. 



9pQa J 



Q 2 



Q 4 



pa 



+ J 2 



pa 



(21) 



where 



TV 



31 1 
240 ~ 60^ 



V24 



115 21 9 7 

J2 H JZ ^4 

48 8 10 



+< — G\ 



+ (—7 G uk. G k\;\v) 



739 2 / 9 3 \ T 133 1 T 19 

1 y + 1 y Ji H Ji H .h Ja 

720 15 y V 8 16 V 48 72 30 

293 3 /55 1 \ 131 145 41 1 

y + 1 y Ji Ji H J3 ^4 

240 10 \24 16 / 16 24 30 



7T ^ 



103 
240 
71 
240 



+ 7^V + 



15 y 



5 
24 



48' 



y)Ji 



59 
48' 



-J 2 



31 

24 



Js 



11 T 2 
— J 4 + -J5 
10 5 



1 1 \ T 61 61 29 2 

-- H y Ji H J2 J3 + —^4 H Js 

8 48 / 48 24 30 15 



29 3 / 1 7 \ 31 23 11 2 

1 y 1 W Ji H J2 J3 H J4 -h 

240 10 \24 16 J 48 24 30 15 



(22) 



There are few checks we can perform to verify our calculation. First, although 
the polarization function in eq.(6) has only two invariant parts, we performed 
the calculation directly leaving the indices /i, v free. This means that we have 
obtained and calculated the 6 different possible tensor structure, given in 
appendix B, separately. Therefore, showing current conservation, g^TI^ = 
q u H flu = is a non-trivial check. Another indirect check is that the final result 
is regular at Q 2 = 0. That is, making a small Q 2 expansion one can show that 
all the Wilson coefficients in eq.(18), eq.(20) and eq.(22) are regular. 
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Eq.(ll), eq.(17), eq.(19) and eq.(21) form the complete OPE up to dimension 6 
operators. Since QCD is renormalizable, there are no other power correction up 
to this dimension and the Wilson coefficients will be an asymptotic expansion 
m 1/logQ 2 . In the vacuum, only scalar operators will contribute. However, in 
medium or when the expectation value is with respect to a nucleon state, the 
tensor operators will also contribute. As an application of our result, we will 
apply our OPE to the QCD sum rule analysis for the J/ip in nuclear matter. 



3. 5 sum of tensor parts 



As we discussed before we will take the nuclear matter to be at rest. Moreover, 
we will use linear density approximation to evaluate the matrix elements, 



<■>'=«• + 2^ (JV| ' W (23) 

where (-) is the vacuum expectation value, p is the nuclear density, m N 
the nucleon mass, and the nucleon state \N) is normalized as (N\N) = 
(27r) 3 2u;Ar£ 3 (0). Then, the OPE from dimension 4 and dimension 6 operators 
in eq.(17), eq.(19) and eq.(21) give the following contribution to II. 



nV) = n 4 , 2 (cu 2 ) + n 6j2 (cu 2 ) + n M (cu 2 
p 



2m 



N 



[C G2 G 2 + (C X X + C Y Y + C Z Z) + C G4 G 4 ] , (24) 



where 



„2 



°G2 = 12(c g )2 [9 - (12 + 2y)J 1 + 9J 2 - 6J 3 ], 
C x = TTTrd^-I- 85 - 2 V + (- 70 + 25 f) J i + 365J 2 - 390J 3 + 252J 4 - 72 J 5 ], 

C Y = ^ 2 3 [25 + 48y + (270 - 45y) Ji - 1185J 2 + 1370J 3 - 408J 4 - 72J 5 ], 
m 2 

°z = olf ,,^ n [-95 - 36j/ + (-130 + 752/) -h + 375J 2 - 190J 3 + 16J 4 + 24J 5 ], 
240(<y^) ci 

4 

°Gi = TT^owI 205 - ( 210 + 33 f) J i + 99 J 2 " 262 J A + 240 J 4 - 72 J 5 ], (25) 

where again in our kinematical limit Q 2 = —uo 2 and the scalar parts of the 
matrix elements in nuclear matter to leading density come from the following 
nucleon expectation values, 



a 1 

(N\-^G%G a av \N) = G 2 {p,p v - -m 2 N g, v ), 
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(N\^Gl x Gl x .^N) = X( P » Pv --^ N g» v ), 
(N\^G* K G a uX . M \N} = Y (p, Pu - -m 2 N g, v ), 

(n\^g; k g: x . m \n) = Z (p,p u - \m\ 9lxv \ 



m 



{N\-^Gl K G a VR . al3 \N) = G 4 p^PvPaPp + —(g^ga/3 + 9^9^/3 + 9^9™) 

1 .2 



48 

+PuPoc9nl3 + PuPpg^a + PaP/39^) 



(26) 



Note that here we chose the nucleoli four momentum to be p = (tun, 0, 0, 0). 
We discuss the magnitudes of these nucleon matrix elements in section 4. 



4 Moment sum rule for J/ip at rest 
4-1 moments 



As an application of our result, we will use our result to refine our previous 
work to calculate the mass shift of the J/ip in nuclear medium[21] using the 
moment sum rule [7,25]. The moments of the polarization function is defined 



^Wo) = ^(-^) n n(g 2 )| Q2 ^. (27) 

where in our kinematics, Q 2 = —uo 2 . Direct evaluation of these moments using 
the OPE gives, up to dimension 4 [25,21], 



M n (0 = ^(0 1 + a n {i)a s + b n (^t + c n (O0c , (28) 



where 



A7T 2 (fG 2 ) 

9b = 



9 (Am 2 c ) 2 
2tt 2 G 2 
3 [Am l c Y 



<Pc = ^TTa — o\o m N Pn (29) 
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and 



= ^ 2> (2n 1 +3)M 1)! ( 4 "' 2 '-"" + «-»«<»■ 5'» + I") 
+ 1)(k + 2)(k + 3) /1 f ,_ 2 2-Fi(n + 2, -j n + |j p) 
2n + 5 1 +4) 2 F 1 (n,i, n +|;p) 

) - k(o - ... 4 " ( : + ^ rt2 ^ + ■ f ' n J h ? < 30 » 

3(2n + 5)(l + 2 2-Pi(n, 2,w+§;p) 

with £ = and p = The factors multiplying the condensate in <\>i and 
4>* in eq.(29) are defined such that at n — > oo c„ ~ 6 n . Moreover, at this limit, 
one notes that the contribution from dimension 4 is proportional to 0^ + 0c> 
which in our kinematical limit originates from the following operator form. 



,4 , ±4 167T 2 (f GIqG%q) 
^ + ^ = — ( 4m 2)2 > ( 31 ) 

where the operator does not have the trace part of the 00 index. That is, 
the leading mass shift is proportional to the color electric field squared E 2 
and the effect coming from the magnetic field squared B 2 disappears. This is 
consistent with the non-relativistic picture at the infinite quark mass limit [16], 
where the Zeeman effect is higher order in a s compared to the Stark effect. 

The dimension 6 operators contribute to the moment as follows, 



+yn(O0' + ^n(O0' + i 74n(O<], (32) 

where 



.3/ 3-1080 (4m;?) 3 
come from the scalar operators and 



3 • 1080 (4m2)3 
2\ 4^(fG^G^} 



(33) 



' 9m 2 N \ 4?r 2 X p 
2 J 3 ■ 1080 (4m 2 ) 3 2m N 

'3m 2 N \ An 2 Y p 
2 ) 3 • 1080 (Am 2 c f 2m N 
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'3m 2 N " 



4tt 2 



94 



(10m 



NJ 



3 • 1080 (4ml) 3 IrriN 
An 2 G 4 p 
3-1080 (Amlf2m N 1 



(34) 



are from tensor operators. Note here again that we have put in the prefactors 
in eq.(33) and eq.(34) such that in the large n limit, all the Wilson coefficients 
become the same. Specifically, the Wilson coefficients are, 



Sn(0 
94n(0 



where 



<7 n (0/n(0, 0; 120, -310, 258, -72) 

s„(0 + <r n (0/„(-120, -60; -300, 1320, -720) 

s„(0 + M0/n(-70, 25; 245, -80, -6) 

Sn{0 + <r„(0/n(270, -45; -1305, 1680, -666) 

Sn(0 + a n (O/„(390, -225; -1245, 880, -306) 

s„(0 + <r n (0/„(-210, -33; -21, 48, -18), 



n(n + 2) 



l + (2n + 5)(2n + 7) 



and 



/n(ci,c 2 ;a 2 ,a 3 ,---,a fe ) 



= ci,c 2 ,a 2 ,a 3 ,---,a fc 



(n + 3)— 



-h P ) 



-2(2n + 7) 2Fl i" +1 ^'"^f 

(n + 3)(n + 4) ^("+Y^ 

(n+3)(n+4)(n+5) 2-Fi("+V 



-§;p) 



2! 



2 Fi(n,i r 



-§;p) 



(n+3)(n+4)-(n+fc+2) 2-Fi("+l, 



(fc-1)! 



2 F 1 (n,i 1 n+|;p) / 



(35) 



(36) 



(37) 



These functions can be read off from the polarizations which are expressed in 
terms of linear combinations of the Jjv(y)'s. See appendix G. 



4-2 estimation of matrix elements 



Here, we summarize the parameters and matrix elements appearing in our sum 
rule. As in ref. [21,25], we will choose the normalization point to be Q 2 = Ami; 
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i.e. £ = 1. Hence, we will use[21] 



« s (8m^) = 0.21, m c = 1.24 GeV. (38) 

For the matrix elements, we will use linear density approximation in eq.(23) 
with m N = 0.93 GeV and the nuclear matter density to be p = 0.17/fm 3 . 



4-2.1 scalar operators 

There are one scalar operator in dimension 4 and two in dimension 6. 

(1) (fG*) p 

The density dependence of the scalar gluon condensate is obtained from 
the trace anomaly relation in the chiral limit which to leading order in 
a s is, 9^ = —^-G a ^ u G a pu . Taking the nucleon expectation value, we find, 

(^G 2 } p =(^G\- 8 - m %p 

~ (0.35 GeV) 4 ^1 - a ^~ S j ■ (39) 

We used m° N = 750 MeV for the nucleon mass in the chiral limit [26] to 
estimate 

(2) G^G' m;m )„ = (-^,){g'f)„ 

Where we have used the equation of motion in eq.(3) to rewrite the 
gluon operator in terms of quark operators. 

This is a four quark operator. Hence to estimate the nuclear matter 
expectation value, we use ground state saturation hypothesis, where one 
can factor out the four quark operator in terms of independent two quark 
operator and take their ground state expectation value[7,27]. This is a 
generalization of vacuum dominance hypothesis(VDH) in the vacuum[7]. 
Hence, 

N 2 -1 r 

(j2) ^ = ~ 2-4^vJ ( ^ Tr[7 ^ ] 

+ (97.9) p(<?7A<?> p Tr [7^7 a 7m ]} 

= (-3) [g^(qq)l + 2 (qi P q) P (qi P q} P - g PP (qi\q) P {qi\q} f 
= {-\)[m? P -mxq) P (qixq) f 
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3 [ m u + m d J 

= -|<??>g (l " a 2 £j ~ -(0.24 GeV) 6 f 1 - a 2 ^ , (41) 



with 



a 2 = ? v ~ 0.69. (42) 

{m u + m d ){qq)o 

In getting this, we first note, 

(?7^)p = (5 = (q i q) p u fl oc p^ , (43) 

where we have introduced the four vector Up to represent the nuclear 
matter, which in our case is u — (1,0,0,0). We used this in eq.(41) and 
neglected terms proportional to p 2 . Other values taken in eq.(41) are, 
(qq) = (-0.23 GeV) 3 and = 0.09/0.014[5]. Finally, we have, 

, a. „„ „„ / 1 



( ^ G^ u G^ a . av ) p ^ (g j ) p 



1 V/Ml-^ 



4tt 2 / V Po 

= (0.23 GeV) 6 ^1- a 2 yj , (44) 

It should be noted that here we have used a larger value of a s = 0.7 
compared to the perturbative value in eq.(38)[10]. 

(3) (^GpuG^u-aa) P 

Using the identities in the appendix, one can rewrite the operator in 
terms of the three gluon operator and the four quark operator. 7^2 (g 3 G 3 — 
g 4 j 2 ) p - For the four quark operator, we use the previous result. For the 
density dependence of the three gluon operator, we assume that the fol- 
lowing ratio is constant to linear order in density. 

This, together with the previous results on the density dependence of the 
gluon condensate and the four quark operator, we find, 



,a x „„ „„ 1 



GG%G^«*)p = ^(9 3 G 3 - g 4 j 2 ) 



7T 



2tt 2 







= (0.60 GeV) 6 f 1 - a 3~\ (46) 
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where 



4-2.2 twist- 2 operators 

For the twist- 2 gluonic operators, we have 

G 2n = -H) 2 "- 2 -A 2 G n , (48) 

7T 

where 

i 

A ( ,;iir) = -2 j<l.r.r- '•(!(.,■. ,r). ( 19) 

o 

We take A%(8m 2 c ) ~ 0.9 and ^(877$ ~ o.02 [28]. 
^.2.3 twist- 4 operators: dimension 6 



Here, there are three independent operators given in eq.(26) and contributing 
to 



(1) (iV|^ A G^JiV> 

This operator can be considered to be the second moment of the gluon 
condensate. For this operator we assume 

(N\^G a KX G: x .^\N) = {p,p u - -g^m 2 N ) x (-)-^(N\-^G 2 \N) 

. 1 2 \ 0.02 / 8 

= \ViiVv - ^9^m N ) x — • \+2m N -m N 

= (P»Pv - ^9^m 2 N ) x (—)m N m° N . (50) 

Therefore, 

X = (^-)m N m%. (51) 

(2) (N\G, K G uX . M \N) = (-^)(g*j^ Ju ) p 

Here, we have used the equation motion again to change this operator 
to four quark operator. However, this operator is traceless and symmet- 
ric. Consider the nuclear matter expectation value of this operator, it 
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becomes, S.T.(g 4 jffi) P = (g 4 jffi) P ~ \g^(g 4 3x3x) P - Assume the ground 
state saturation hypothesis again, 

TV 2 — 1 

(3;3i) P = - 2 ; A2 . N {(qq)l^hM 

+{qq) P {qi.q) P ^[i»i K iv + i K i P i»] 

+ {qinq) P {qixq) P ^[ 1 K llll x lu \} 

= (-g) [gpv(qq) 2 P + ^{qi P q)Mivq) P - g pv (qi\q) P (qi\q) P 

(3 a X) P = (-\) [m? P - 2(qixq) P (qixq} P ] ■ (52) 



So 



2 _ 1 

S-T-(3^fu) P = -^{qi P q) P (qivq) P + -g P Ml\q) P (qi\q) P - (53) 



Using eq.(43), 



^.r.o^^^^tiV- 4 "^) (54) 



Hence, 

y = (55) 

(3) (N\^G« K G: X . M \N) 

Here we use the equation of motion once. Then one finds that the 
operator has similar structure to a twist-4 quark gluon mixed operator 
appearing in the twist-4 contribution in deep inelastic scattering [2]. 

(N\G° K G a KX . M \N) = {N\gj K [D v , G^N) 

~ (N\iq{D„ *F l/X } 1 x l5 q\N) 

= ^{N\mw 5 {D v , *F^ K }u\N) 
1 3 

= {v»Vv - -^g^m 2 N ) ■ -K 9 U . (56) 

Here, K 9 is a parameter introduced in ref. [3,4] to fit the available deep 
inelastic scattering data. The fit gives, -0.3 GeV 2 < K g u < -0.2 GeV 2 . 
With this one finds, 

Z ~ (0.097 GeV) 2 . (57) 

In Table 1, we summarize values of the re-scaled condensates in Eq. (29), (33) 
and (34) in the vacuum and at normal nuclear matter density. 
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Table 1 

Expectation values of normalized gluon operators in the vacuum and at normal 
nuclear matter density 



Element 


<t>i 


4>i 


<t>t 


<t>t 


in vacuum 


1.67 x 1(T 3 





1.66 x 1(T 7 


1.65 x 1(T 8 


in matter 


1.58 x 10~ 3 


-1.53 x 10~ 5 


1.28 x 1(T 7 


6.50 x 10~ 9 


change 


-9.70 x 1(T 5 


-1.53 x 1(T 5 


-3.88 x 1(T 8 


-1.00 x 10~ 8 


Element 










in vacuum 














in matter 


1.67 x 1(T 8 





4.27 x 1(T 10 


-1.69 x 10~ 9 


change 


1.67 x 1(T 8 





4.27 x 1(T 10 


-1.69 x 10~ 9 



4-3 numerical analysis 



The polarization function IT satisfies the following energy dispersion relation. 



9n 1 f , ImII(s) ,_„ s 
n (u 2 ) = - / ds K -±. 58 

TV J S — UJ 

In the vacuum, the spectral density (ImIT(s)) consists of the J/ip pole, the 
contribution from the excited states ip' and the DD continuum. The contribu- 
tion from the J/ip and the low-lying resonances can be approximated by delta 
functions inside the dispersion integral. This is so because the J/ip mass lies 
below the continuum threshold and is dominated by electromagnetic decays. 
This is also true in nuclear matter for a J/ip at rest. The inelastic channels 
opening due to the scattering of the J/ip + N — > A c (2.28) + 5(1.87) is for- 
bidden by kinematics. All other processes are OZI rule violating. Hence, the 
delta function approximation for the J/ip is also good in the nuclear medium. 



Imn(s) = fS(s - m 2 m ) + f'S(s - m%,) + .. (59) 

The second term just represents a generic excited state contribution. Substi- 
tuting eq.(59) into the dispersion relation, taking the moments and taking the 
ratio between the (n — l)th and the nth moment, one finds, 



Mn^Q 2 ) 2 ^ (1 + 6, 



™ 2 JI* + Q 2 ) >< ( ) ( 60 ) 



where 
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m\, + Q 



^7k^ • (61) 



Substituting the vacuum values for the excited states, one finds that the ratio 
1 t S 'l~ 1 goes to 1 for n > 5. Hence, the mass is determined from the relation in 
eq.(60) at Q 2 = 4m* 



_ M n ^{Aml) 2 



(62) 



As before[21], we analyze this equation as a function of n. In Fig. 2(a) , we plot 
the previous result, which includes the contribution only up to dimension 4[21]. 
In Fig.2(b), we plot the present result which includes the total dimension 6 
contribution. As can be seen from the comparison, the minimum occurs again 
at similar n value and the change from the vacuum results are also similar. 
We avoid fine tuning of the bare charm quark mass m c to fit the vacuum J/ ip 
mass to its vacuum value, because we are only interested in the shift of the 
J/ip mass, which is almost independent of this fine tuning. 





Fig. 2. The mass of J/ijj in GeV determined from the plot of eq.(62) for different 
n at £ = 1. We compare the result in normal nuclear matter (solid line) with the 
vacuum result (dashed line), (a) refers to the previous calculation, which includes 
only dimension 4 operators, (b) is the present calculation which includes all the 
dimension 6 operators. 

Comparing the two graphs, one notes that the minimum occurs at the same n 
value and the graphs looks similar. Comparing the changes from the vacuum 
curve and the medium curve at the minimum point, we find 



Am 



j/1> 



-4 MeV. 



(63) 



1 The determination at different value of Q 2 was found to be not significant [22] 
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This mass shift is smaller than our previous calculation including dimension 4 
only. The main reason for a smaller mass shift compared to including dimen- 
sion 4 only is as follows. In the vacuum, the dimension 6 operators tend to 
cancel the dimension 4 operators [10]. This tendency is not only true but more 
effective in the medium. Therefore, including dimension 6 effects in medium 
would effectively correspond to a smaller change in the dimension 4 operators 
in our previous analysis in ref[21]. This implies a smaller mass shift. Below, 
we will try to elaborate on this point and to explain the interplay of each 
contribution. 

• In Fig. 3, we plot the contributions from dimension 4 and dimension 6 op- 
erators to the moments both in the vacuum and in nuclear medium. As 
expected, the dimension 6 contributions are smaller compared to dimension 
4[21]. Nevertheless, one notes that both in the vacuum and in the nuclear 
medium, their contributions are opposite to each other so that the contri- 
bution from dimension 6 operators tend to cancel the contributions from 
dimension 4 operators. 





Fig. 3. Contributions from dimension 4 and dimension 6 operators to the moments 
both in the vacuum and in the nuclear medium in GeV~ 2n . 

• The mass shift is coming from the changes in the contributions from each di- 
mension. Hence, for comparison, in Fig.4, we plot the contributions coming 
from each dimension in the vacuum and in the nuclear medium. One notes 
that effectively, each contributions becomes smaller in nuclear medium. 
Numerically, the ratio between the changes in dimension 4 operators and 
dimension 6 operators are roughly given by (dimension 4 : dimension 6 
-3:1). 

• The changes in dimension 4 operators are dominated by the changes in the 
scalar gluon condensate [21]. This is evident by comparing <f)f and in table 
1. The comparison between <pf and 0^ is enough to compare its contribution 
to the moments, because the Wilson coefficients multiplying them are the 
same in the large n limit. Moreover, b n = c n oc n 3 and b n — c n oc n such that 
b n ~ c„ at moderate n. 
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That is not quite so for the dimension 6 operators. Here, the Wilson 
coefficients multiplying 0^,..'s given in eq.(34) are also defined such that 
they are equal in the large n limit and are proportional to n 5 . However, the 
difference in the Wilson coefficients differ by only one power of n, that is, 
the difference between Wilson coefficients are of order n 4 . Therefore, the 
Wilson coefficients are substantially different at the value of our interest, 
which is smaller than n = 10, and become similar only at very large value of 
n > 100. Hence, we will analyze each contribution from dimension 6 more 
in detail to identify the important contributions. 

In Fig. 5 (a), we have plotted the density dependent part of the scalar, 
twist-2 and twist-4 contributions. As can be seen, the important contribu- 
tions are the scalar operators. And among the scalar operators, the contri- 
butions coming from <f>® is more important than that from (ffi s . It should be 
noted that this is so because of the relatively large Wilson coefficients at 
small n values. The relative importance among the scalar operators are the 
same also in the vacuum. The next important set of operators are the twist- 
4 operators. In Fig. 5 (b), we plot the contributions from twist-4 operators, 
which has no vacuum expectation values. The most important contribution 
is coming from (G a G a . a/3 ) , which contributes to 0^. Here, the Wilson co- 
efficients are similar at least among the twist-4 operators and the relative 
importance can be read off from Table 1. 
• We finally compare the dominant contributions from dimension 4 and di- 
mension 6 operators to the moments. In Fig. 6 (a), we show the total con- 
tributions of the scalar dimension 4 and scalar dimension 6 operators to 
the moments. As can be seen from the figure, the dominant contributions 
from dimension 6 operators have opposite sign from dimension 4 operator 
and reduces its contribution. This reduction is enhanced for the density de- 
pendent part. This is shown in Fig. 6 (b). The reason why the change in 
the dimension 6 operator are greater is clear. The gluon condensate eq.(39), 
which is the dominant dim 4 operator, changes in nuclear matter only by 6 
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Fig. 5. (a) shows the contributions from the density dependent part of the dimension 
6 scalar , twist-4 and twist-2 operators to the moments, (b) shows the individual 
contributions of the twist-4 operators. 



0.1 



-0.1 



-0.3 



-0.5 





-&--ET"~1j] 


i — -g M_o — a — a — a — a — a — «- 

c e G 2 ^\ 




* fl G, a G,p 








.2 

A A j 




(a) 





Q. 



Q 



0.03 



■P 0.02 



CJ 0.01 



-0.02 



a- a G 3 G 3 ,, 

(iv nv:aa 



(b) 



Fig. 6. (a) shows the important contributions to the moments coming from scalar 
operators of dimension 4 and dimension 6. (b) shows the density dependent part of 
(a), i.e. the density dependent changes in the condensates. 

%, whereas, the scalar four quark condensate eq.(41), which is the dominant 
dim 6 operator, changes by almost 70%. Therefore, although the Wilson co- 
efficients are smaller for dimension 6 operator, the density dependent part 
of dimension 6 operator are large such that it cancels the the contribution 
from dimension 4 operator non-trivially. 

Hence, the reason why the mass shift gets smaller compared to just tak- 
ing into account dimension 4 operators is because the density dependent 
changes coming from dimension 4 and dimension 6 operators tend to cancel 
each other. Therefore, taking into account dimension 6 operators would ef- 
fectively be equal to a smaller change in the dimension 4 condensate in the 
previous result [21], which would have given a smaller mass shift. 
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5 Conclusion 



In this paper we have calculated the OPE of the correlation function between 
two vector currents made of heavy quarks up to dimension 6 operators with 
any tensor structure. The formidable task was to categorize and calculate the 
corresponding Wilson coefficients of dimension 6 twist-4 gluon operators. This 
is a first attempt to establish the three independent twist-4 gluon operators 
in terms of operator basis. 

Using this result, we have applied our OPE to analyze the mass shift of J/ip 
in nuclear medium using QCD moment sum rules for the heavy quark system. 
This is a generalization of our previous result [21], where we calculated the 
mass shift using the OPE only up to dimension 4 operators. Unfortunately, the 
nucleon expectation values of the dimension 6 operators are not as reliable as 
the dimension 4 operators. Nevertheless, using an order of magnitude estimate 
for the matrix elements, we find that the mass of J/ip would decrease by 
about 4 MeV in the nuclear medium. This is 3 MeV smaller than the previous 
result on including only dimension 4 operators, and shows that the dimension 
6 effect is about 40% correction of the dimension 4 effects and goes in the 
opposite direction. This result seems consistent with the notion that the higher 
dimensional correction in the vacuum QCD sum rule for the heavy quark 
system goes like (G 2 /mf) K in the r n = ^"^""^ with alternating signs with 
if [29]. This also seems to be true in medium and the true mass shift is expected 
to lie between —4 and —7 MeV. 

The resulting value of mass shift, is also consistent with the more recent esti- 
mates using a totally different approach [17, 18]. This is a result for a J ftp at 
rest with respect to the nuclear medium. However, since we have calculated 
the OPE for a general external four momentum, our results can be easily and 
reliably generalized to study the moving <//■?/> [24] and also to finite tempera- 
ture[30], which would be also interesting in relation to the ongoing discussion 
of J /if) suppression in RHIC due to a comover model. 
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A Identities for spin-2 dimension 6 gluon operators 



The spin-2 dimension 6 gluon operators in eq.(4) are not independent. The 
following relations holds among them. 



^ kX;/j, uk;X 



fiK\X k\;u 
{iK\v kA;A 



( /-to, I ria \ _ /~ia ^L/^" 1 f^ a 



/-ia (~ia 
fin\vX kX 



2 

-{ U uk;Xv + 9 J ^uk^vX^kX 



■ {-)gfG^ v + G a ^ K . x G a KX . u 
Or, 

/k;A sA;i/ yiK\v kX;X 

2 kA;(i"kA;i' ' ^ yiKiu^ kA;A 

kX\Xv [ik - - [^kX^X + 9 J ^kX^Xu) ^ fx 

' kX,v [ik,X y^fiu fiK;v kX;X 

r~ia r~ia ( r~ia i „ eabcr~ib r~ic \ r~ia 

fj,K;Xn uX - - { U »k;kX + 9 J ^^Xk) ^ v 



vX 



9G 3 



yw 



[ik;X vk;X 



yiK\K^ J uX\X ' 2 kX;[i kX;v ' yin\v kA;A 

: ~G\\ K -\\G a UK = —g ^2/ a6c G^ A G'^ A + j°. K — j^.^j G a ul , 
■. {-2)90% + g 2 fX - 9f K G a ^ 

^ kX;^ kX;u ZLT t «t;i/ Ur KA;A 5 J/Jz/ ^k^kAjA 

✓-fa ✓-la , 2 -a -a , 

fxn;v kX,X y -V-^ ^kAj^^kA^ 



(A.l) 
(A.2) 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



Using these, one can reduce the operators in eq.(4) to the three independent 
operators in eq.(5). 



B Current Conservation in the Polarization Operators 



Let us consider the polarization of vector currents defined in eq.(6). When the 
current is Hermitian, i.e. = j^, the contribution from spin-2 dimension 4 
operator can in general be written in terms of the following four independent 
terms. 



IV (g) = ^ 



i , gVg/^ 2 , 

+ Q4 "a? + -Q2 



yiu 



no 3 



Q 



Q 2 



4 

a/3 



(B.l) 
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Imposing current conservation, on finds, 



f) 2 =f) 1 4- ft 4 



(B.2) 



Similar relations holds for higher dimensional operator. This is the operator 
form of showing the existence of two independent polarization directions in 
medium for eq.(6). 

Summarizing similar relations for different spins and dimensions, we have 
(1) dimension 4 and spin 2 



n*'„ 2 (g) = 



■ L liv ' r < 1 \1p1ll- L pv T QpQv^pfj,) + 9pv 2 i7 pCT + 4 V-pa+Jpa) 



(2) dimension 6 and spin 2 



(<2 2 ) 2 



V + T^ilplP-hv + Ipl^pp) + dpv—^-Jpa + 



(Ipa + J P-) 



(3) dimension 6 and spin 4 



(Q 2 )3 



r 4 j 1 r 4 i t-4 \ , ^p^o" r 4 

1 K \p,v + T^WP^'sAp^ + 1p1v l K\pp) + dpi' o2 J K Ap<r 



~l~ q4 V nXpa ' J K\pa) 



C Propagators 



Here we summarize the quark propagator in the presence of the external gauge 
field in the fixed point gauge [23]. 



iS(p) = J d 4 xe ipx iS(x,0) 

= iS m (p)+ / d 4 xe ipx g I d 4 ziS m (x- z)i A(z)iS {0 Hz) 



+ ••• 



(C.l) 



and 
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iS(p) = J d A xe~ ipx iS(0, x). (C.2) 
In the fixed point gauge, we write the field in terms of covariant operators, 



+^yX a xpx p {D a DpG Pfl (0)) + --- (C.3) 



Collecting terms, we can write the full propagator in terms of gauge covariant 
fields. A few symbols are used for convenience: 



' _ • ® <- — ■ 9 (n a\ 

and 



{a, (3} = 7 Q - 7/3- 



m f> — m Tb — m 



{a, (3, a} = 7 a - 7^- 7^- 



$ — m $ — m $ — m $ — m 
• 1 

and P(a, (3, 7...) means sum of all possible permutations in a,/3,j.... 
Here we simply list the propagators. 
(1) iSW(p)=iS<®(p) 



(C.5) 



— m 

(2) iS G (p) = iS G (p) 

= iSW(j>) Unp ■ \x a G aP ) iS^{p) = \G a p{^ a} 
= iSW(p) i( 7/3 • \ x a G a/3 ) iS(°)(p) = -^G a/3 {a, /?} 

(3) iScp(p) = iS G 2(p) 

= iSW(p) i( lf3 ■ h a G aP ) iS^ip) i( 7<T ■ \x p G pa ) iS^{p) 
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; ~~^G a {3 G p a a, a, p} + {(3, a, P(a, p)} 
iS^(p) i( 7/3 ■ l - x a G a/3 ) iS(°)(p) z( 7ct • \ x p G pa ) iS<®(p) 
-^G a pG pcr [{a, (3, p, a} + {P(a, p), (3, a}J 



(4) As for G 3 part, 
• iS G s(p) 



{13, a, <t, p, A, k} + {/3, a, a, A, P(p, k)} 
+ {/3, <t, P(a, p), A, k} + {/3, <t, p, A, P(a, k)} 
+ {P, <r, a, A, P(p, k)} + {/3, a, A, P(a, p, k)} 



iS G 3{p) 



iS (0) (p) H-yplxaGcp) iS<°>(p) i( la \x p G pa ) iS<°>(p) i{~ix\x K G KX ) iS<°)(p) 



-G a aG pa -G K \ 



{a, f3, p, a, k, A} + {P{a, p), f3, a, k, A} 
+ {a, 13, P(k, p), a, A} + {P(a, k, p), (3, a, A} 
+ {P(a, p),f3, k, a, A} + {P(a, k), (3, p, a, A} 



(5) As for DG part, 
• iS DG {p) 



■-iS (0 Hp) ■ \x a x p D a Gp a ) iS^ip) 



-D a G Pa {a,P(a,(3)} 



iS DG {p) 



= iSW(p) i{ la ■ \x a xpD a G Pa ) iS(°)(p) 
= l -D a G l3a {P(a,p),a} 



(6) As for D 2 G part, 
• iS D 2 G {p) 



= iS {0 \p) ■ \x a xpx v D a D p G va ) iS<°)(p) 

o 

= --D a D p G VIT {j,P(a,P,v)} 
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iS D 2 G {p) 



■iS (0) (p) i( 1<T ■ \x a x p x v D a D p G va ) iS(°)(p) 



= +-D a D f3 G l/a {P{a, a} 
(7) As for (DG) 2 part, 



-D a Gf3 K DpGcr\ 



{K,P(a,f3),\,P(p,a)} 



+ {«, a, A, P(/3, p, cr)} + {k, /3, A, P(a, p, a)} + {k, A, P(a, (3, p, a)} 



iS(DG) 2 (p) 



{P(a,(3),K,P(p,a),\} 



+ {P(a, 0, p), k, a, A} + {P(a, /3, cr), k, p, A} + {P(a, /3, p, cr), k, A} 



(8) As for GD 2 G part, 

• iS GD 2 G (p) 

= iSW(p) i( 7/3 • ^x a G a/3 ) iS(°)(p) *( 7ct • ixpf K x AJ D pJ D K G A(T ) 
= ^G a pD p D K Gxa [{/?, a, <r, P(p, k, A)} + {/?, cj, P(a, p, k, A)} 

• iS GD 2 G {p) 

= iS ( -°\p) i( 7/3 ■ l -x a G aP ) iS^(p) i(la ■ \x p x K xxD p D K G Xa ) iS<® 
1 G a pD p D K Gxa (3, P(p, k, A), a} + {P(a, p),0, P(k, A), a} 
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+ {P(a, k), /?, P(p, A), a} + {P(a, A),/?, P(k, p), a} 
+ {-P(a, «, A), /3, p, a} + {P(a, p, A),/3, k, cr} 
+ {P(a, k, p), A, cr} + {P(a, p, k, A), /3, cr} 



(9) As for D 2 GG part, 
• iS D 2 GG (p) 
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D p D K G Xa G a p [{a, P(p, «, A), /?, a} + {a, P(«, A), /?, P(a, p)} 
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+ {a, p, (5, P(k, A, a)} + {a, (5, P(p, k, A, a)} 
+ {a, P(k, p),(3, P(a, A)} + {a, «, /?, P(a, A, p)} 

+ {fi, P(A, p), /?, P(a, «)} + {a, A, /?, P(p, a, «)} 

• iS D 2 GG (p) 

= iS^(p) i( 7tr • ix p x K x A D pJ D K G ACT ) iS^fr) t( 7/ j • ^«G a/3 ) iS(°)(p) 
= r^D p D K G Xa G af3 [{P(p, k, A), a, a, /?} + {P(p, /e, A, a), a, /?}] 



D spin structure of operators of dimension 6 



In computing the Wilson coefficients we used the following reduction of the 
Lorentz indices to the spin-2 operators. 



fabcsia fib fic _ A — An — An -I- A r 

J ^ uu af3 pa ■ rL pa'^vpf3a ■ rL fip'^upaa ■ rL i/a'^ppf3a "T ■ rl u/3'^fj 

-^fip^uaa/3 -^pa^uap/3 -^-i/p^paa/3 -^-ua^pap/3 

H~ AapC/3p<jv ^aaCpppu A-fip^apav ^fia^appv (^•1; 



^ /n^l H2V2;af) ~ ^ afi^ pxp2VW2 

~^~-P/3piCap2ViV2 Pfiv\Cetp2P\V2 
Jap2^Ppi V2V\ Jau2^fSpip J 2i^i 

a pi C f3 fi2 V1V2 ^ T av\C'fip2P\V2 

^i- 1 p\P2^fiv\av2 i- , u\p2^fip\av2 

Lp 1 i/2df3uiap2 L 'v\V2^ '/5 'p\a.p2 

~^~^-p\P2^ficiv\V2 ^-v\p2^ , a.fip\V2 

p\p2^ ftV20LV\ r ^v\p2^'fiv20ip\ 
~ ^ p\V2^fip2Civ\ r ^v\V2^'fip2Cip\i (D-2) 



where 



Cafipv — QaflQpu QavQpfi) 

daf3pv = 9a(39pu ( D . 3 ) 

and, if we take 3 operators, ddggl = G® X G® X .^ V , ddgg2 = —g 2 j^j®, and 
dc%3 = G* K G a KX . Xv as our basis, we get 
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A = {ddggl + 2dc%3)/4, 

K = 2P = 2 J = (13ddggl + ddgg2 + 19ddgg3) /80, 
Q = W = (-27ddggl + ddgg2 - 61ddgg3)/160, 
L = (Addggl + 3ddgg2 + 7ddgg3)/A0, 
M = (29ddggl + \3ddgg2 + 47ddgg3) /160, 

T = (-Gddggl + 3c%</2 - 13dc%3)/40. (D.4) 



As for the spin 4 part, we have a simpler reduction: 



G \xivfi \i2V2\af) ^aj3fiifj,29u!i/2 ^ ' <xfSv\V29 Hi (J.2 

(D.5) 



where 



- 2 G ViK G V2K;ap- ( D -6) 



E Integrations with respect to Feynman Parameter 

In general, after Feynman integration, one can write the polarizations in terms 
of linear sums of J's defined in eq.(13). 

n(g) = Q' a [a J + aiJi + a 2 J 2 + ■ ■ ■ + a k J k ] . (E.l) 



To reduce the polarization function into this final form, we use the following 
steps and identities. 

After the Feynman integral, the polarization function will be a sum of Iran's: 



I% n (Q 2 , ™ 2 ) = I dx X) w = IT{Q\ m 2 ) (E.2) 

1 [m 2 + Q 2 :r(l - x)] N 



o 



We then follow the following steps. 

(1) I mn can be expressed in terms of 1% = I™ using the following identity. 

X n + (1 - X) n = 1 - | nClx"- 1 ^ - X) +„ C 2 ^- 2 (l - z) 2 + • • • +n C* fc X™- fc (l - x) k 

+ ---+ n C n - 1 x(l-x) n - 1 } 
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= 1 -„Ci{x"- 1 (l - x) + x(l - x)™- 1 } 
- n C 2 {x n -' 2 (l - x) 2 + x 2 (l - x) n -' 2 } 

= 1 - n d X(l - X)|x™- 2 + (1 - XT' 2 } 
-nC 2 X 2 (l - x) 2 {x n - 4 + (1 - X)"- 4 } 

(2) Then we can reduce 1% to 1% = ijv using 7$ = ^ (^/v-i _ m 2 !^ 1 ) = 



1 rn- 1 1 jn—1 

Q2 1 N-1 y 1 N-l 



(3) We then introduce the dimensionless function Jn(u) = In, where 
y = Q 2 /m 2 . 

(4) Finally, we use the recurrence relation. 



2 4N — 6 
y J N (y) = — — - + — — - Jjv-i(y) - 4Jjv(y) 



iV-1 JV-1 
(5) We can also integrate explicitly, 



(E.3) 



Jn(v) 



(27V - 3)!! 



(AT-1)! 

where s — 1 + 4/y. 



/s-lX" , \ - (fc-1)! fs-lV 



(E.4) 



F An example: evaluation of a diagram 



Diagram la: We show the evaluation of the Feynman diagrams here. The 
first diagram in Fig. 1 is taken as an example. 



d 4 k 

d 4 k 



Tr[7 M S< >(fc + 9)7„5 G3 (fc)] 



Tr 



7m" 



1 



G h pcr G% x t a t b t c 



{a, (3, p, a, re, A} + {P(a, p), f3, a, k, A} 



/ (2tt) j ' AH (/- <>! 

+ {a, /3, P(re, p), a, A} + {P(a, re, p), /3, <r, A} + {P(a, p), 0, re, <t, A} + {P(a, re), /3, p, <r, A} 

1 



1 * / rate r i) \ / d 4 k 



-Tr 



g 4 V" cp „c Ki a. i j ^4 ( k + g y 2 _ m 2 

X {a, (3, p, a, re, A} + {P(a, p), /3, tr, re, A} + {a, f3, P(re, p), a, A} 
+ {P(u, re, p), 13, a, A} + {P(a, p), /3, re, <t, A} + {P(a, re), /3, p, er, A} 

""h A v pC^ iOLG p A va C^ a pfi + A a pCp^ i( j u ^-aa-C/^^ipi/ J ^/3p c Q/[(Ji' ~^ ^^(r^a/ipi/ j 

y ( 27r )' 



4 (fc + g) 2 -™ 2 (fc 2 -m 2 ) 7 



7m (#+ tf+mfrv 
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: j(Jf + m) 7a (# + m) 7/3 (/(S + m) lp {}i + m)-y a (jt + m) TK (# + m) 7A (# + ; 

/ !_ fl fW 

7 (2tt) 4 (k + q) 2 -m 2 (k 2 -m 2 ) 7Jla( ' 

' r(8) 



(2^) 4 



,.(> 



dx 



[{k + (1 - x)q} 2 - S 2]8 r(7 ) 



/!«(*) 



1 

i f , f d 4 l 1 r(8) f 28 o /116 2 220 2 2 \ lfi 

o L 

+ (-^m 4 + ^-m 2 q 2 t 2 - ^qH^ I* + (20m 6 - 60mVt 2 + 60mV< 4 - 20< ? 6 t 6 ) « 2 }a m! , 

+ j(-640t 3 + 896t 4 )i 4 + (640m 2 i 3 - 768m 2 t 4 )Z 2 

-128m 4 (t 3 - t 4 ) + 256mV(t 5 - i 6 ) - 128<? 4 (t 7 - t 8 ) ^q^q a q0A a/3 
f 140 , fi /464 2 400 2 2 80 , i\ ,4 

+(-148m 4 t - 112m 4 t 2 + 96m 2 q 2 t 3 + 296m 2 g 2 * 4 + 116g 4 t 5 - 184<? 4 t 6 )/ 2 

+(40t + 24t 2 )m 6 - (56 9 2 * 3 - 72<jr 2 t 4 )m 4 + (-8g 4 t 5 + 72<jr 4 t 6 )m 2 + 24q 6 (t 7 - t s )^g^q a q fj A al3 
f / 14 o\ , fi , 464 2 560 2 2^4 

+| ( — t- 4or J r + (-— m * + — m 2 r)/ 4 

+(148m 4 t - 176m 4 t 2 - 4Wm 2 q 2 t 3 + 472m 2 <? 2 i 4 + 268<? 4 t 5 - 296<? 4 t 6 )Z 2 

-40m 6 (t - t 2 ) + 120mV (t 3 - t 4 ) - 120mV(t 5 - * 6 ) + 40g 6 (i 7 - t 8 )} (c? M g a A al , + g„ ?a A aM ) 

7T 2 / 5 52 \ 

+9t**>Qaq0A al3 — y- + y- (2 + j/)Ji - 9J 2 + yJ3 - 8J 4 J 

/ gg gg 

+qnq v q a qpA a p ^ y - (8 + 4j/) Ji - 36J 2 + y J3 - 8J4 



7T 2 / 

+{q tl q a A aI/ + qvq a )A alJ . — y 



7r / 23 ,2 1 w 41 20 

3 --»-(5 + 5»W- T A + 1 -A 



) 



(F.l) 



G Moments in terms of Hypergeometric Functions 2 -Pi 



In order to obtain moments written in terms of hypergeometric functions 2 -Pi 
from polarizations eq.(E.l), we expand the polarizations in terms of Q 2 /m 2 . 
Here we divide ao and a\ into ao = d\ + d^y and ai = c\ + C2y. The other 
higher coefficients, 02, 03, • • • , a& don't contain y, where y = Q 2 /m 2 . 



n(q) = Q- 



ao + (ci + C2J/) 



[l + x(l-x)Q 2 /m 2 ] 



[l + x(l-a;)Q 2 /m 2 ] 2 
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+a k 



[1 + x(l - x)Q 2 /m 2 ] k 



Q- a 



o-o 



f <2 2j 

(ci + c 2 y) I dx2_^ {-x(l - X )Y —2j 
n J=° 



i 

/oo 
{-x(l - 

3=0 



i 

/oo 

3=0 



(fc-l+j)! Q 2j ' 



(k - m?i 



=0-5: 



3=1 



i i 

n / dx(-l)V(l -z) J ^V + a 2 / d:r(-l)V'(l-x) j (l + j)^- + 
m l 3 I m z 3 



/ <M- U'.r'iJ - r)'^ — -C, 2 j dx{-\)'x' ' ( 1 - .()•' 



/ 



1 



: {k-l+j)\ Q 2j 

(k - m?3 



dx {x(l — x)} J 



(_Q2 )3 -a/2 



,2j 



j=a/2 



+ ^(l + j)(2 + j)- 



c(l-x)|ci+o 2 (l+j) 



(k - 1)! 



ft — 1 

nc+i)} 



+ C2 



(G.l) 



Now the n's moment M„ can be obtained from above by differentiation. 



(-) 



a/2_ 



( 4m 2)a/2+n (! + £)n+l 



if r(j)r(n+f +2) 1 .a + 3 ? 

°2 — tt^n 2-Fl(" + 1, — 



+ 



a 3 r (5) r («+f + 3 



S+3\ 2F 1 {n + l,--,n 



2' ' 2 1 + 5 
3 a + 3 £ 



2 '1 + 5 



2! T(n+^) 
a4 r(|)r(n+f +4) 5 a + 3 £ 

"3T r(n + ^) 2Fl(n + 1 '-i' n+ ^'— } 



+ • 



+ (*-!)! r(n+^) ^(» + 1 '- fc +2' B+ — + 



i (-) a / 2 r(i/2) 



( 4m 2)a/2+n ( X + £)n+l 2r(n+ £+3) 



j^£(l/2) / a w \ 
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/ (n+§)/0! 2 Fi(n + l,i,n+^±2;p) \ 

-2(2n + a + l) 2 Fi (rt + 1, ± , n + 2±1 ; p) 
(n+f)(n+2±2)/l! 2 Fi(n + 1, -i, n + 2±2 ; p ) 
X (n+f)(n + ^i±2)(n+S±4)/2! 2 Fi(n + 1, -f ,n + 2±2;p) 

\ (n + §)(n + 2±2) • ■ ■ {n + § + k - l)/(fc - 1)! 2 Fi(n + 1, -fc + § ,n + 2±2; p) ) 

(G.2) 

where 



' 1 T[l/2] r[n+(a+2)/2] p ( , 1 1 , a+3 . g ^\ 
2 T[n+(a+3)/2] 2^ H" t 1, 2 , H t 2 , 1+ £ ) 

As an example, when £ = and a = 6, we have, 





-( 







(G.3) 



M n (Qg,g = 0): 



(-1) 



(4m?)" 



i r(i)r(n + 4) r(l)r(n + 3)(n+|) 

C l- — ^\ +C 2 (-2) ; — h 



2 ri 



if r(i)r(» + 5) Q3 r(i)r(n + 6) 
"U a2 r(„ + f) + 2- r(n+f) 



(-1) 2 n+3 (n + 2) 



( 4m 2)n+3 (2n + 7) 



2)! ( \ 
— I ci,c 2 ,a 2 ,a 3 ,- • • ,a fe I 



M° 
(4ml) 



2 5 ?r 2 -n(n + 2) 
3 (2n + 5)(2n + 7) 



( ci,c 2 ,a 2 , • • • , a k ^ 



r(»+f) 

g fc r(i)r(n + fc + 3) 
"(fc-i)i r(n + f) 

/ (n + 3)/0! \ 

-2(2n + 7) 
(ra + 3)(n + 4)/l! 
(n + 3)(n + 4)(n + 5)/2! 



. (n+2+k}\ , v I 

\ (n + 2)! /( k - 1 y- ) 

( (n + 3)/0! \ 

-2(2n + 7) 
(n + 3)(n + 4)/l! 
(rt + 3)(n + 4)(n + 5)/2! 



(G.4) 



where Mj = 3x 2 n (rt + l)(n - l)!/(47r 2 (2n + 3)!!(4m 2 )") = A„(£ = 0). 



H Consistency check 



As discussed in the text, there are few checks we can perform to confirm our 
calculation. The first is the current conservation, which we checked explicitly. 
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The second is the regularity at Q 2 = 0. This can be checked from the following 
equation. 



n(,) = 



i 

oo „ 

E (-) a/2_1 / 

,■ in J 



dx {x(l — x)} 3 



.! (-Q 2 y-«/ 2 



i2j 



j=a/2 



e(l-x){ 



ci + a 2 (l +j) 



+ ^ ( i +J )(2 +j ) + - + ^n^^} + - 

i=i 

holds only after making the following checks: 
• j = 



(H.l) 



= d\ + ci + a 2 + a 3 H + a k 



3 = 1 



c 2 + d 2 + J dx {-x(l - x)} {ci + 2a 2 + 3a 3 H + /ca fc } 



J = 2 



i i 
= J dx {-x{l - x)} c 2 + J dx {-x(l - x)} 2 |ci + ■ ■ ■ + fc(fc 2 +1) a fc | 



J = 3 



1 1 



fc(k + l)(fc + 2) 



-ife 



± 

= /" dx{-x(l - x)}^ 2 c 2 



+ y dx{-x(l-x)}2 1 |ci + --- + i ^— ^ 2 - J -a k j 

So we know the summation begins at j = a/2. 

Substituting the values for the coefficients multiplying J's, we have ex- 
plicitly checked that the above constraints are satisfied in our calculation. 
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